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Abstract. This paper is devoted to the thermally activated dynamics of the 
capillary condensation. We present a simple model which enables us to identify the 
critical nucleus involved in the transition mechanism. This simple model is then 
applied to calculate the nucleation barrier from which we can obtain informations 
on the nucleation time. We present a simple estimation of the nucleation barrier 
in slab geometry both in the two dimensional case and in the three dimensional 
case. We extend the model in the case of rough surfaces which is closer to the 
experimental case and allows comparison with experimental datas. 
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1. Introduction 

When two surfaces are brought together in a condensable vapor near saturation, a first 
order phase transition from gas to liquid occurs at small gap width provided that the 
solid wets the solid substrate i.e. has a contact angle smaller than 90°. Macroscopic 
considerations predict that the condensation occurs for distance between the solid 
surfaces H less than a critical distance He verifying 

Ap A/i 2{jsv - 1sl)/Hc (1) 

where Ap = pi — pg is difference between the bulk densities of the liquid and the gas 
and A/i = psat — /i is the (positive) undersaturation in chemical potential, with psat 
the chemical potential at bulk coexistence [1] . If the vapor is assumed to be a perfect 
gas, then Ap « ksT \n{Psat / Pvap) = ksT \n{l / RH), where RH is the the so-called 
relative humidity. At standard ambient conditions for water {jlv = 72 mJ.m~^, 
PL ~ 3.10^^ m~^, H = 40%), we obtain He ^ 2 nm. Capillary condensation is 
usually invoked to interpret adsorption isotherms of gases in mesoporous media [2]. 
This transition is now well understood and documented, both from the experimental 
[3, 4, 5] and theoretical point of view [6, 7]. 

On the other hand, the problem of the dynamics of the transition has receive 
very little attention. Experimentally, only indirect informations on the dynamics are 
available in the litterature. Experimental studies of capillary condensation using the 
Surface Force Apparatus (SFA) technique systematically show an important hysteresis 
in the interaction force between two substrates when the separation of the surfaces 
is first decreased and then increased. This large hysteresis points out the strong 
met ast ability of the gas phase when H < He, which persists over macroscopic times. 
Recent experiments with SFA have studied the growth of the liquid meniscus after 
the nucleation [8] but no attention has been given to the nucleation time which is 
much larger than this growing time. Experiments measuring the cohesion inside 
divided materials may provide indirect information on the dynamics of the transition 
too [9, 10, 11]. Theoretically, lattice-gas simulations showed that the (topologically 
equivalent) drying transition occurs via the creation of "tubes" connecting the two 
wetting films [12]. 

Beyond these results, a theory proposing a mechanism for the dynamics of the 
capillary condensation is still needed. 

In the following, we show how to construct the critical nucleus for capillary 
condensation. First, a simplified model in the slab geometry keeping only the 
main ingredients for capillary condensation will be considered. The latter has both 
advantages to allow tractable calculations and to capture the essential features of the 
involved physics. In a second part, we show how the natural roughness of the surfaces 
on a nanometric scale can be taken into account on the dynamics of the transition. 
Applications to the adsorption kinetics in a granular medium shall be discussed. 

2. The slab geometry 

In a first step, we restrict our attention to a system confined between two perfectly 
smooth and fiat solid surfaces, and in contact with a reservoir of temperature T and 
chemical potential p. 

Let us consider the situation in which planar liquid films of varying thickness e 
(e < H/2) develop on both solid surfaces. Following Evans et al. [13, 6], the grand 
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potential of the system may be written 

n = -pvVv - PlVl + 2jslA + 2-fLvA (2) 

where Vy (resp. Vl) is the volume of the gas (resp. liquid) phase and A is the surface 
area. Using Vl = 2Ae, Vy = A{}i — 2e) and pv — Pl — ApA/i, one gets 

Au;{e) = ^{n-n{e = 0)) = ApA/i 2e (3) 

Note that in the complete wetting situation ft{e = 0) can be identified with 
the grand-potential of the system filled with the gas phase only. The situation 
e = H/2 corresponds to the opposite case where the two liquid films merge to fill 
the pore. The grand potential thus exhibits a discontinuity in e = H/2 corresponding 
to the disappearance of the two liquid- vapor interfaces, and its value is reduced by 
2'y LY A. When e = H/2, expression (3) must then be replaced by Auj{e = H/2) = 
—ApAjii{Hc — H), where He is the critical distance defined in eq. (1). One may 
note that the minimum of the grand potential corresponds to a complete filling of the 
pore by the liquid phase when < i^c, as expected. If we now allow deformation 
of the interfaces, i.e. the thickness e is now a function of the lateral coordinates, the 
corresponding cost has to be added to the grand potential. We assume also a mirror 
symetry of the interfaces, so that one finds in this case: 

^^tot = Jlv^Alv + J dSAu;{e) (4) 

with AQtot = ^({e}) — and AAly = Aly — A is the excess L — V area. The 
integration in the last term runs over the solid surface. 



2.1. The 2D case 

Let us consider first the 2D case. Within the small slope assumption, AAlv — 
J dx 7Ly|Vep, extremalization of the grand potential leads to the following Euler- 
Lagrange equation for e(x), where x denotes the lateral coordinate: 

d^e dAuj(e) 

We look for solutions satisfying e = and de/dx = at infinity. We can choose 
e{x = 0) = H/2 to fix the origin. The complete solution, depicted in fig l.a, can 
be obtained in the form of parabolic branches with a spatial extension Xc = \/HRc 
where Rc = He/ 2. Let us note that the cusp in the solution in x = stems from 
the discontinuity of Auo in e = H/2 resulting from the assumption of an infinitesimely 
narrow liquid- vapor interface. Condensation thus occurs through the exit at ion of 
short wavelength fiuctuations, in agreement with the simulations results for the drying 
transition [12]. The corresponding energy of the nucleus (per unit length in the 
perpendicular direction) can be calculated by integration of eq. (4) : 

Afit = l(A;xAp7i^)i/2H3/2 (6) 

It is easy to check that AVi^ corresponds to a saddle-point of the grand-potential. 
It is greater than both free energies of the gas and liquid phases. Moreover AQ/^ is 
smaller than the free energy of any other configuration maximizing the grand potential 
since it is the only solution of finite extension. We just point out that the parabolic 
solution obtained above is the small slope approximation to the circle with radius of 
curvature Rn. 
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Figure 1. a : Picture of the critical nucleus for capillary condensation in two 
dimensions and perfect wetting case (0 = 0). The radius of curvature of the 
meniscus is equal to Rc = He/ 2, which is only approximatively verified within 
the small slope assumption. See text for details, b: Picture of the critical nucleus 
in three dimensions and perfect wetting case. See text for details. 



We mention that the prediction for AQ^ in eq. 6 is in agreement with numerical 
simulations results, using a Landau- Ginzburg model for the grand potential of the 
system together with a non conserved Langevin dynamics. Full details of these 
simulations are given elsewhere [14]. 

These results can be generalized to the partial wetting case. The only difference 
is that the contact angle on the surfaces is now non vanishing and obeys Young's 
law, jsv — ISL = iLV cos6. To leading order on H/Hc, one gets [14] : 

A(^t ^ 27Ly sin (9i^ (7) 



2.2. The 3D case 

The previous approach can now be directly generalized to the 3D case. The physical 
understanding of the results can be however greatly simplified if one realizes that 
maximization of the grand potential, eq. (4), leads to two mechanical equilibrium 
conditions : the usual Laplace equation, relating the local curvature to the pressure 
drop ^LVf^ = — A/iAp; and the Young's law which fixes the contact angle of the 
meniscus on the solid substrate according to jlv cosO = jsv — Jsl {i-e. = in the 
perfect wetting case). These non-linear equations cannot be solved analytically in 3D, 
but one can easily understand that the corresponding critical nucleus takes the form 
of a liquid bridge of finite lateral extension i?*, connecting the two solid surfaces (see 
fig lb). This finite extension results physically from the balance between a "surface" 
contribution Al^i ^ {ApA/j.H—2{jsv—1sl))7tR'^ which drives capillary condensation, 
and a linear contribution Aft2 ^ 27rjLvRH specific to the 3D case which tends to 
close the bridge. Maximization of the free energy gives a finite extension i?*, yielding 
for the free energy barrier : 

2i^sv-7sL)He-H ^ ' 

Full details in the 3D case shall be given in a forthcoming paper [15]. 



3. The rough case 



Although a lot can be learned from the perfectly fiat slab geometry, the latter is 
certainly too idealized to account for the kinetics of adsorption in "real" experimental 
systems. In particular, very slow logarithmic depedence have been measured on 
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various static properties of granular media in the presence of humidity (see fig 3 and 
text hereafter) [9, 10]. As we shah show hereafter, these logarithmic time dependence 
may be understood by taking into account the influence of roughness on the dynamics 
of capillary condensation. Let us consider a simple model consisting of two surfaces 
facing each other and rough at the nanometric scale, as depicted on fig. (2. a). As 
emphasized in the introduction capillary condensation typically occurs in pores of 
nanometric size. We thus have to consider the roughness of the surfaces at the 
nanometer level. Without loss of generality, one may consider that one of the walls 



Vd 




a b 

Figure 2. a: Typical representation of two rough surfaces. Note that we consider 
the roughness at the nanometric scale. 

b : Schematic representation of an asperity, is the excess volume of the defect, 
ttd the area of the defect. 



is perfectly flat. When roughness is present, there is a broad range of gaps between 
the surfaces. In particular, there are regions where the two surfaces are in close 
contact. In such regions, condensation should take place on a very short time-scale. 
Thus at "early times", one has to consider a set of wetted islands, which we shall 
consider as independent. Once these islands have formed, they should grow up to a 
point where the distance between the surfaces is equal to i^c, so that a meniscus of 
radius Rc = He / 2 cos forms at the liquid- vapor interface, allowing for mechanical 
equilibrium. 

In doing so however, the wetted area has to overcome unfavorable regions where 
the distance between the two surfaces is larger then He. Let us consider a specific 
jump over such a "defect", as idealized in fig. (2.b). We denote by Cd the "averaged" 
gap inside the defect {cd > He), and by ad its area. The free energy cost for the liquid 
bridge to overcome this defect is approximatively given by 

An^ ad (A/iAp Cd - 2A7) 

= VdA/j^Ap (9) 

where Vd is the excess volume of the defect, Vd = ad (e^ — He). We can thus estimate 
the time to overcome the defect as 

One may expect the defects to exhibit a broad distribution of excess volume 
Vd, so that the activation times r are accordingly widely distributed. After a time 
only the defects with an activation time r smaller than t have been overcomed. 
Using eq. 9 and 10 these have an excess volume Vd which verifies Vd < Vdmaxif) = 
kBT{AiiAp)~^ ln(t/ro). At a time t, the number of overcomed defects is then typically 
N(t) = Vdmax{t) / vq where vq is the typical width of the distribution of excess volume 
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Figure 3. a : Evolution of the mass m of a pile of glass beads of radius smaller 
than 50 fim as a funtion of the logarithm of the resting time t in hours. Note 
that the time is comprised between a few minutes and 2 weeks. The temperature 
is fixed at 31 ± 0.1 °C. The relative humidity is fixed at 68% by the salt method 
described in reference [10]. The straight line is the best linear fit of the datas. b 
: Evolution of the tangente of the maximum stability angle Om of a an assembly 
of glass beads as a function of the logarithm of the "resting" time t in seconds 
(divided by the cosine of this angle which comes from geometrical arguments). 
This angle is measured in a cylinder. The full experimental setup is described in 
[10]. (• ) : RH = 3%, (■) : RH = 43%. The straight fines are the best finear 
fits of the datas. 



of the defects. Now, once a liquid bridge has bypassed a defect, it fihs locahy the 
volume surrounding the nucleating site and the wetted area increases by some typical 
(roughness dependent) amount SAq. The time dependent wetted area can thus be 
written : 

Similar expressions with logarithmic dependence on time, can be found on other 
quantities, like the time dependent adsorbed amount, or the adhesion force between 
rough surfaces. 

These logarithmic dependence have been observed in two kinds of experiments. 
In the first one (see fig. 3. a), we have measured the evolution of the mass of a 
sample of glass beads with a diameter smaller than 50 micrometers at fixed humidity 
{RH = 68%) as a function of the resting time t [16]. The glass beads where first 
dried at high temperature and put at the fixed humidity controlled by the saturated- 
salt method described in [10]. The evolution of the mass fits well with a logarithmic 
behavior, as described on eq. 11. On the other hand, the cohesion force resulting from 
condensation of liquid bridges has been probed in a granular medium by measuring 
the maximum angle of stability as a function of resting time. As shown on figure 3.b, 
the latter exhibits a slow logarithmic dependence in agreement with eq. 11 [9, 10]. 

The authors would like to thank J. Crassous, J.C. Geminard and H. Gayvallet 
for many interesting discussions. This work has been partly supported by the PSIVIN 
at ENS-Lyon, the IVIENRT under contract 98B0316 and the franco-british program 
ALLIANCE (contract 99041) 
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